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Chapter  1 


INTRODUCTION 

Standard  Base  Supply  System 

The  Standard  Base  Supply  System  (SBSS)  Is  an  auto¬ 
mated  Inventory  accounting  and  control  system  designed  to 
provide  timely  support  to  base  level  activities.  The  system 
uses  the  UNIVAC  1050-11  computer  for  storage  and  maintenance 
of  records  and  for  the  generation  of  management  reports. 

The  system  as  a  whole  consists  of  both  manual  components  and 
interfacing  computer  programs.  The  system  consists  of  four 
major  functional  processes--item  accounting,  accounting  and 
finance,  file  maintenance,  and  management  reporting  (24«l-2). 
When  using  the  UNIVAC  1050-11  on-line  computer  system,  all 
records  affected  by  a  given  transaction  are  updated  at  the 
time  of  transaction  input.  The  system  provides  remote 
devices  in  work  centers  on  base  to  insure  timely  communica¬ 
tions  between  man  and  machine  (24il-5). 

The  SBSS  was  designed  to  help  speed  the  flow  of 
materials  from  many  sources  to  the  base  user  organizations. 
The  SBSS  inventory  is  stocked  by  material  from  Air  Force 
depots,  the  General  Services  Administration  (GSA) ,  Defense 
Supply  Agency  (DSA)  and  through  local  procurement  (19il). 

The  SBSS  is  driven  by  user  demand  actions. 
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The  S3SS  automatically  provides  outputs  to  other 
management  systems.  Those  outputs,  such  as  requisitions, 
go  to  wholesale  suppliers  such  as  Air  Force  Logistics  Com¬ 
mand  (AFLC)  and  the  Defense  Logistics  Agency  (DLA) .  The 
SBSS  imposes,  Air  Force  wide,  standard  organization,  pro¬ 
grams  and  procedures.  The  organization  is  designed  for  ease 
of  customer  support  and  efficiency  of  internal  operation 
(23«l-6).  The  entire  SBSS  is  an  extension  of  the  basic 
supply  requirements  to  order,  receive,  store,  and  issue 
property. 


Inventory 

Inventories  are  used  throughout  the  world  in  many 
organizations,  both  civilian  and  military.  Inventories  can 
be  classified  in  many  ways.  The  following  classification 
is  one  convenient  wayi 

1.  Production  Inventories i  raw  materials,  parts 
and  components  which  enter  the  organization’s  product  or 
service  in  the  production/ transformation  processes. 

These  may  be  either  special  items  manufactured  to  user 
specifications  or  standard  items  purchased  "off  the 
shelf." 

2.  MRO  inventories!  maintenance,  repair  and  oper¬ 
ating  supplies  consumed  in  the  production/transformation 
processes  but  which  do  not  become  part  of  the  product  or 
service. 

3.  In-process  inventories!  semifinished  products 
found  at  various  stages  in  the  production  operation. 

4.  Finished -goods  inventories!  completed  products 
ready  for  shipment  and/or  use  [l4il89J. 

Regardless  of  their  classification,  inventories 
serve  the  primary  purpose  of  decoupling  successive  stages 
in  the  production-distrlbution-consuraption  chain  (I3il). 
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Inventories  allow  production  decisions  and/or  supplier 
demand  decisions  to  be  made  independent  of  supplier  procure¬ 
ment  decisions  where  demand  is  independent  or  near  independ¬ 
ent. 

Inventory  Control 

Inventory  control  is  a  vital  element  in  the  manage¬ 
ment  of  system  resources.  Development  of  analytical  tech¬ 
niques  and  computer  capability  have  combined  to  transform 
inventory  control  into  a  critical  function  requiring  profes¬ 
sional  managerial  skills  (14il87). 

The  significance  of  the  inventory  control  function 
is  clearly  demonstrated  by  the  level  of  resources  involved. 
Studies  have  shown  that  a  commercial  firm’s  inventory  com¬ 
monly  constitutes  anywhere  from  15  to  25  percent  of  its  total 
invested  capital  (14il87),  In  its  everyday  operations,  in¬ 
ventory  control  is  often  a  computerized  operation  within  a 
carefully  defined  and  controlled  structural  framework 
(14 i 22 5).  Sound  management  in  the  area  of  inventory  control 
is  vital.  Management  must  try  to  balance  the  various  risks 
of  low  inventory  with  its  associated  risk  of  stock  out, 
production  halts,  back  orders j  and  high  Inventory  with  its 
risk  of  high  carrying  cost  and  increased  risk  of  obsoles¬ 
cence  (23il42).  To  achieve  this  management  objective  a 
forecast  of  future  demand  must  be  made.  This  forecast  is 
based  on  uncertain  economic  cycles,  customer  demands,  and 
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advances  In  technology.  In  light  of  this  uncertain  environ¬ 
ment,  inventory  control  must  answer  how  much  to  order  and 
when  to  order  it  (l3i2). 

Generally  speaking,  three  basic  systems  have  been 
employed  in  controlling  inventories i  (1)  the  cyclical  order¬ 
ing  system,  (2)  fixed  order  quantity  system,  and  (3)  mate¬ 
rial  requirements  planning  system  (14il94).  The  first  sys¬ 
tem  is  a  time-based  system  involving  scheduled  periodic 
reviews  of  the  stock  level  of  all  inventory  items.  The  sec¬ 
ond  type--the  fixed  order  quantity  system--is  based  on  the 
order  quantity  factor  rather  than  on  a  time  factor.  The 
material  requirements  planning  system  concept  provides  a 
way  of  looking  at  the  management  of  production  inventories 
in  a  dependent  demand  environment  (14 i 194),  Regardless  of 
the  system  used,  management  must  have  an  accurate  assessment 
of  what  will  happen  in  the  future  in  order  to  make  its  deci¬ 
sions  (13 i 2). 

Forecasting 

In  almost  any  activity,  having  perfect  knowledge  of 
the  future  would  be  an  immense  help.  Rather  than  the  cer¬ 
tainty  of  perfect  knowledge,  the  inventory  manager  finds 
himself  in  an  uncertain  environment.  He  can  do  little  to 
change  that  environment  so  he  seeks  management  techniques 
which  will  aid  in  reducing  the  levels  of  uncertainty.  All 
forecasting  techniques  have  in  common  the  basic  goal  of 
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reducing  the  uncertainty  in  one's  expectation  of  the  future 
(4il).  As  Robert  L.  Sims  says,  "It  should  be  obvious  that  I 

to  forecast  is  to  be  wrong.  Only  in  the  most  fortuitous 
cases,  and  these  seldom  occur,  can  a  forecaster  guess  the 
future  exactly  £2l«l].M  Today,  the  economics  of  spiraling 
prices  and  austere  budgets  make  reliable  estimates,  or  fore¬ 
casts,  especially  critical  ( 3 1 22) •  Forecasting  for  and 
acquiring  spare  parts  is  an  important  facet  of  the  material 
support  system.  Forecasts  are  used  to  develop  dependable 
systems  and  to  formulate  decisions  concerning  the  techno-  j 

logical  and  economic  feasibility  of  retaining  a  system  in 

j 

use  (Hi  15).  The  exponential  acceleration  of  spares  cost 
provides  a  compelling  reason  for  basing  requirements  corapu- 

j 

tation  upon  the  best  available  forecasting  methodology 
(3.22). 

Standard  Base  Suoolv  Forecasting  System 

The  Air  Force  Logistics  Command  (AFLC),  throughout 
its  Air  Logistics  Centers  (ALCs)  or  depots,  buys  and  dis¬ 
tributes  centrally  procured  items  of  supply  used  on  mis¬ 
siles,  aircraft,  and  other  equipment.  Each  depot  buys  and 
stocks  items  in  specific  federal  supply  classes  (FSCs)  and 
serves  as  the  primary  source  of  supply  for  Air  Force  bases 
and  other  activities  ( 24 « 11-1), 

The  SBSS  is  driven  by  user  demand  actions.  When  a 
base  needs  an  item  a  requisition  is  submitted  to  SBSS.  If 
the  item  needed  is  not  currently  in  stock,  a  request  for 
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that  item  is  placed  on  order.  The  request  may  be  filled 
through  the  normal  operating  cycle.  In  some  instances,  spe¬ 
cial  orders  may  be  needed.  Orders  for  resupply  of  the  SBSS 
are  logged  as  due-ins  and  are  maintained  until  the  material 
is  received  (19il). 

Repair  parts  and  supplies  with  expendability, 
repairability,  recoverability  category  (ERRC)  designator 
XB3  will  be  stocked,  using  a  variable  economic  order  quan¬ 
tity  (EOQ)  stockage  concept.  An  XB3  item  is  an  expendable 
item,  i.e.  an  item  consumed  in  use  or  incorporated  into 
another  assembly.  The  EOQ  model  is  applicable  when  the 
quantity  of  items  ordered  arrives  in  the  inventory  at  one 
point  in  time  and  when  the  demand  for  the  item  has  a  con¬ 
stant,  or  nearly  constant,  ratej  and  the  cost  of  the  item 
is  relatively  low  (2i497).  The  SBSS  requirements  determina¬ 
tion  may  be  divided  into  the  range  model  and  a  depth  model. 
The  range  model  refers  to  the  procedure  used  to  determine 
if  an  item  is  to  be  stocked  at  the  base  level.  The  depth 
model  is  used  if  the  item  is  to  be  stocked  to  determine  how 
much  to  order  and  when  to  order  it.  The  depth  model  is  based 
on  the  classic  EOQ  formula! 

EOQ  »  V  2DA/IP  (1.1) 

where,  D  *  annual  demand  rate 

A  ■  cost  per  order  (currently  used  figure,  $4,54) 

I  ■  annual  inventory  carrying  rate  (currently  used 
figure,  26  percent) 

P  ■  item  unit  price  ( 19 1 2 ) 


This  formula  balances  the  cost  of  ordering  with  the  cost  of 
holding  inventory.  The  derived  order  quantity  will  minimize 
these  total  variable  costs  (19j2). 

The  Air  Force  application  of  the  EOQ  concept  pro¬ 
duces  a  variable  requisition  objective  by  determining  the 
EOQ  in  consideration  of  the  number  of  demands,  daily  demand 
rates,  and  a  stockage  priority  code.  The  purpose  of  the 
variable  EOQ  stockage  concept  is  to  prevent  the  premature 
stockage  of  items  based  on  erratic  or  indefinite  demand  pat¬ 
terns  (24 j 11-2),  The  basic  EOQ  model  has  been  enhanced  to 
better  meet  the  variable  stockage  objectives  (VSO)  of  the 
SBSSj 

EOQVSO  -  V  2*DDR*VS0»A/(I.P)  (1.2) 

where,  A  =  cost  per  order  ($4.54) 

I  =  annual  inventory  carrying  rate  (26  percent) 

P  =  item  unit  price 
DDR  =  daily  demand  rate 

„  VSO  =  number  of  days  in  the  EOQ  computation  (19i3) 
Another  important  concept  incorporated  into  the 
SBSS  is  that  base  supply  levels  are  managed  on  a  reorder 
point  basis.  That  is,  each  time  the  stock  position  of  an 
item  reaches  or  falls  below  the  established  reorder  level, 
some  type  of  supply  action  must  be  taken  to  bring  stocks 
up  to  quantities  authorized  (24 i 11-2),  The  reorder  point 
is  a  combination  of  the  order  and  shipping  time  quantity 
(OSTQ)  and  the  safety  level  quantity  (SLQ) .  The  OSTQ  is 
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that  quantity  required  by  a  base  to  permit  uninterrupted 
replacement  from  the  external  supply  source.  The  OSTQ  is 
given  by i 


OSTQ  =  DDR* OS T 


(1*3) 


where,  DDR  =  daily  demand  rate 

OST  =  average  order  and  shiptime  in  days,  based  on 
the  item  source  and  priority  ( 19 *  3) 

The  safety  level  quantity  is  a  variable  quantity,  in  units, 

computed  to  provide  protection  from  stockouts  during  the 

reordering  process.  The  SLQ  is  given  by« 


SLQ  =  C  V  3* OSTQ  (1.4) 

where,  C  =  the  safety  factor  (typically  set  to  1,  which 
implies  an  84  percent  service  effectiveness, 
given  the  assumption  of  normally  distributed 
demands) 

3  =  lead  time  demand  variance/mean  ratio 
OSTQ  =  order  and  shiptirae  quantity  (19«3) 

The  reorder  point  (RP)  is  given  byi 


RP  =  OSTQ  +  SLQ  C 19 i 3)  (1.5) 


The  maximum  desired  inventory  position  is  called  the  requi¬ 
sition  objective  (RO) i 


RO  =  INT  (EOQvso  +  OSTQ  +  SLQ  +  0.999)  (1.6) 


where,  INT  ■  integer  result  of  the  following  computation 
EOQvso  “  enhanced  EOQ  model 

OSTQ  =*  order  and  shiptime  quantity 
SLQ  ■  safety  level  quantity 
0.999  ■  factor  added  to  the  EOO  level  to  adjust  to 
the  next  highest  unit  (19*4) 
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It  can  be  seen  that  each  of  the  three  variables  in  the 
requisition  objective  computation  are  a  function  of  the 
daily  demand  rate.  The  daily  demand  rate  is  the  forecast 
measurement  for  the  SBSS.  If  inaccurate  estimates  for  DDR 
are  made,  errors  may  result  (19.3) . 

'./hen  an  item  is  ordered  the  daily  demand  rate  is 
revised.  The  SBSS  maintains  status  on  the  cumulative  recur¬ 
ring  demand  (CRD)  and  the  date  of  first  demand  (DOFD).  Each 
time  an  item  is  ordered,  the  number  of  units  ordered  is 
added  to  the  CRD.  The  revised  DDR  is  them 

DDR  =  CRD/max  (180,  Current  Date  -  DOFD)  (1.7) 

The  minimum  of  180  days  usage  is  assumed  so  as  not  to  over¬ 
stock  items  that  have  been  recently  added  to  the  stockage 
list  (I9i3).  The  number  of  demands  for  each  ERRC  code  XB3 
item  will  be  recorded  in  six-month  increments  up  to  a  total 
of  three  increments  or  18  months;  i.e.,  the  current  and  two 
past  six-month  increments.  When  18  months’  demands  have 
been  accumulated,  the  oldest  six-month  increment  will  be 
dropped  and  a  new  six  months'  accumulation  will  begin.  The 
criteria  for  establishing  demand  levels  for  EOQ  items  (ERRC 
code  XB3)  will  vary  dependent  upon  the  stockage  priority 
code,  but  will  not  be  less  than  three  demands  per  year 
(24.11-3). 
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Problem  Statement 


The  literature  search  has  shown  that  a  significant 
error  variance  exists  in  the  forecasting  method  employed  in 
the  Standard  Base  Supply  System  when  compared  to  the  actual 
demand  experienced.  The  current  forecasting  system  used  by 
the  SBSS  is  a  combination  of  exponential  smoothing  and  mov¬ 
ing  averages.  The  literature  suggests  that  this  technique 
may  not  provide  the  best  estimate  of  future  demand, 

A  critical  measurement  used  in  SBSS  is  forecasted 
demand.  A  forecasting  approach  which  reduces  forecast  error 
has  potential  impact  beyond  its  effect  upon  economic  order 
quantity  (EOQ)  determination.  Accurate  demand  data  in  the 
context  of  the  SBSS  affect  average  inventory  levels  and 
holding  costs,  and  Improve  order  and  shipping  time  (lead- 
time),  thereby  impacting  safety  stock  levels  and  premium 
transportation  requirements  for  priority  back  orders. 

This  thesis  explores  whether  a  multiple  model  fore¬ 
casting  technique  can  improve  demand  forecasts. 

Research  Objectives 

The  objectives  of  this  thesis  arei 

1.  Evaluate  the  multiple  model  forecasting  tech¬ 
nique  as  an  alternative  to  the  Standard  Base  Supply  System 
forecasting  technique, 

2,  Determine  which  forecasting  approach  will  pro¬ 
vide  the  most  accurate  estimate  of  future  usage. 


3.  Determine  the  effectiveness  of  implementing  the 
proposed  system. 

4.  Recommend  actions  based  upon  the  results 
obtained . 

Research  Justification 

This  thesis  applies  the  research  of  Mitchell  and 
Garland  at  the  wholesale  level  of  demand  to  retail  level 
forecasting.  Those  authors  used  simulated  and  limited  actual 
data  to  conclude  that  the  multiple  model  forecasting  tech¬ 
niques  actually  forecasted  more  accurately  than  the  single 
method  currently  used  in  the  Air  Force's  wholesale  level 
D062  system  (13«42).  This  thesis  had  the  sponsorship  of  the 
Air  Force  Logistics  Management  Center  (AFLMC) ,  Gunter  Air 
Force  Station,  in  a  continued  effort  to  improve  present 
inventory  control  methods. 

Scope 

This  thesis  will  utilize  a  research  data  base  main¬ 
tained  by  AFLMC  and  containing  two  and  a  half  years  of  data 
from  Dover  Air  Force  Base,  Delaware.  The  data  base  contains 
approximately  15,000  line  items  (National  Stock  Numbers). 

The  data  were  tested  using  several  different  forecasting 
techniques.  The  results  were  compared  to  results  obtained 
by  previous  research  on  the  SBSS. 
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Chapter  2 


LITERATURE  REVIEW 


Background 

Major  Richard  Lombardi  of  the  Stocking  Policies 
Division,  Air  Force  Logistics  Management  Center,  summarized 
the  current  dilemma  of  logisticians  faced  with  forecasting 
requirements  to  meet  mission  demands » 

The  Standard  Base  Supply  System  (SBSS)  demand  fore¬ 
casting  is  an  unorthodox  prediction  scheme  which  may 
lead  to  suboptimal  forecasting.  This  is  attributed  to 
the  manner  in  which  the  daily  demand  rate  (DDR)  is  up¬ 
dated  by  six  month  adjustments  in  the  SBSS  cumulated 
recurring  demand  (CRD;  and  the  date  of  first  demand 
(DOFD)  fields.  The  net  effect  of  those  adjustments  is 
to  convert  the  forecasting  to  a  form  of  exponential 
smoothing  in  which  the  value  of  the  smoothing  constant 
will  vary  depending  upon  date  of  demands.  This  varia¬ 
tion  may,  depending  upon  when  demands  are  experienced, 
provide  more  weight  to  past  demands  at  the  expense  of 
current  demands,  thus  dampening  the  influence  of  cur¬ 
rent  requirements  [15], 


Multiple  Model  Forecasting  Method 

In  order  to  choose  an  effective  forecasting  system, 
the  forecaster  must  decide  upon  a  model  which  is  most  appro¬ 
priate  given  the  conditions  that  exist  at  the  time  of  the 
forecast.  A  variety  of  forecasting  models  are  available 
from  which  the  forecaster  can  choose.  The  selection  of  a 
forecasting  method  depends  on  many  factors i  the  context  of 
the  forecast,  i.e.,  the  nature  of  the  decision  environment 


12 


at  the  time  of  the  forecast,  the  relevance  and  availability 
of  historical  data,  the  desired  accuracy  of  the  forecast, 
and  the  time  period  to  be  covered  (13«ll).  Due  to  these 
factors,  certain  types  of  forecasting  techniques  are  better 
suited  to  a  particular  demand  pattern  than  others.  Since 
in  practice  we  are  seldom  faced  with  a  pure  "classical" 
demand  pattern,  the  forecaster  should  consider  several 
methods  for  forecasting  a  single  demand  pattern  before 
selecting  the  model  best  suited  to  the  actual  demand  pat¬ 
tern. 

The  multiple  model  forecasting  method  used  in  this 
thesis  research  incorporates  forecasting  techniques  recom¬ 
mended  by  a  variety  of  experts.  These  techniques  are  moving 
average,  single  exponential  smoothing,  double  exponential 
smoothing,  and  adaptive  exponential  smoothing.  The  computer 
was  programmed  with  these  simple  forecasting  strategies. 
Through  the  process  of  simulation,  the  computer  will  select 
the  one  best  strateg}  to  forecast  an  item  at  a  given  moment 
in  time.  Whichever  strategy  best  projected  the  most  recent 
completed  period  is  the  one  the  computer  uses  to  project 
the  future  (22 » 3). 

Focus  Forecasting 

A  new  approach  to  inventory  control  using  the  con¬ 
cept  of  declining  computer  costs  per  computation  is  advo¬ 
cated  by  Bernard  T.  Smith  in  his  book,  Focus  Forecastingt 
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Computer  Techniques  for  Inventory  Control  (22).  Focus  fore¬ 
casting  as  implemented  by  Smith  used  a  series  of  simple 
forecasting  approaches.  It  uses  the  powerful  computation 
simulation  ability  of  the  computer  to  pick  the  one  forecast¬ 
ing  strategy  that  will  work  best  for  the  one  item  for  the 
next  period  (22«xii).  Focus  forecasting  goes  through  all 
the  computations  every  time  it  forecasts  a  demand. 

Several  parallel  research  efforts  have  been  accom¬ 
plished  using  various  forecasting  techniques  at  the  wholesale 
level  on  the  AFLC  D062  inventory  control  system  for  expend¬ 
able  items  and  the  AFLC  D041  reparable  asset  management  sys¬ 
tem.  A  summary  of  these  theses  follows. 

Garland  and  Mitchell 

The  purpose  of  their  study  was  to  determine  if  a 
multiple  model  forecasting  technique  could  forecast  demand 
more  accurately  than  the  model  currently  used  in  the  Air 
Force  Logistics  Command  D062  system  for  expendable  items. 
Simulated  and  actual  data  were  used  to  check  the  results. 

The  methods  utilized  in  the  multiple  model  technique  were 
an  eight-quarter  moving  average,  a  four-quarter  moving  aver¬ 
age,  exponential  smoothing,  adaptive  smoothing,  a  least 
squares  fit  and  a  ratio  of  change  between  years  method. 
Results  were  compared  in  terms  of  mean  absolute  deviation 
adjusted  to  show  percentage  change  in  accuracy  compared  to 
the  D062.  The  statistical  test  used  for  comparison  was  the 
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t-test  for  matched  pairs.  This  test  indicated  approximately 
a  17  percent  improvement  in  accuracy  using  either  simulated 
or  real  historical  data. 

Christensen  and  Schroeder 

This  research  effort  compared  the  D041  Single  Moving 
Average  forecasting  method  used  to  forecast  reparable  gen¬ 
erations  of  recoverable  items  with  the  Box  and  Jenkins*  time 
series  analysis  forecasting  methods.  Five  artificially  gen¬ 
erated  stochastic  processes  were  used  to  model  the  possible 
reparable  generations  observed  in  practicei  (1)  a  Poisson 
process  with  a  constant  mean,  (2)  a  Poisson  process  with  a 
decreasing  mean,  (3)  a  Poisson  process  with  an  alternating 
linear  mean,  and  (4)  a  process  whose  values  are  the  sine 
function  of  the  output  of  a  Poisson  process.  The  research 
concluded  that  the  D041  forecasting  method  made  unbiased 
forecasts  for  the  Poisson  process  with  a  constant  mean  and 
the  sine  function,  but  made  biased  forecasts  for  the  other 
three  processes.  Time  series  analysis  forecasting  methods 
were  only  used  to  make  forecasts  for  the  processes  that 
were  found  to  be  biased  using  the  D041  forecasting  method. 
Time  series  analysis  forecasting  methods  made  unbiased  fore¬ 
casts  for  the  processes  whose  means  were  linearly  increasing, 
linearly  decreasing,  and  alternating  linearly. 
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Brantley  and  Loreman 

This  research  effort  examined  the  forecasting  tech¬ 
nique  used  in  the  Air  Force  Logistics  Command  D041  reparable 
asset  management  system.  It  was  hypothesized  that  the  mean 
of  the  absolute  values  of  the  D041  forecast  error  in  prac¬ 
tice  was  equal  to  zero.  This  hypothesis  could  not  be 
rejected.  It  was  further  hypothesized  that  another  time 
series  forecasting  technique  (exponential  smoothing)  would 
also  yield  an  error  distribution  with  a  mean  of  zero.  This 
hypothesis  could  not  be  rejected  eitherj  moreover,  the  vari¬ 
ances  for  the  exponential  smoothing  forecast  error  distribu¬ 
tions  were  less  than  the  D041  forecast  error  distributions 
for  all  four  lead  times  examined  (one,  two,  three  and  four 
quarters) . 


Fischer  and  Gibson 

This  thesis  effort  examined  the  application  of 
exponential  smoothing  to  forecast  demand  for  economic  order 
quantity  items.  The  research  team  stated  that  management  of 
EOQ  items  required  the  use  of  a  demand  forecasting  technique 
to  estimate  future  demand  for  the  purpose  of  establishing 
stock  levels.  The  study  compared  the  effectiveness  of  four 
forecasting  models  that  could  be  used  at  base  level.  The 
moving  averages  method  and  single,  double  and  triple  expo¬ 
nential  smoothing  were  evaluated  using  22  months  of  demand 
history  for  a  random  sample  of  34  EOQ  items  stocked  at  a 


base  consolidated  supply  activity.  Four  statistical  error 
measures  were  used  to  compare  the  accuracy  of  the  forecasts 
generated  by  the  models  for  the  items.  The  study  concludes 
that  the  methods  were  not  significantly  different  in  fore¬ 
casting  EOQ  items  in  terms  of  accuracy,  stability  and  bias. 
The  study  team  recommended  a  reporting  system  be  considered 
chat  would  permit  the  study  of  the  application  simultaneously 
of  various  forecasting  techniques  to  the  management  of  an 
item  and  the  factors  affecting  item  demand  patterns. 

Patterson 

The  purpose  of  this  research  was  to  look  into  alter¬ 
nate  approaches  to  forecasting  demand  for  expendable  items 
in  the  SBSS.  The  forecasting  models  studied  were  single, 
double  and  adaptive  exponential  smoothing.  Sample  items 
were  selected  from  the  actual  demand  records  from  Dover  AFB. 
An  analysis  of  the  various  models  was  performed  by  means  of 
a  computer  program  written  for  each  model.  The  forecasting 
models  were  compared  on  the  basis  of  forecast  error  as  meas¬ 
ured  by  the  mean  absolute  deviation  (MAD) ,  The  forecast 
error  for  the  model  currently  used  by  the  SBSS  was  also 
measured.  The  research  study  concluded  that  the  single 
exponential  smoothing  method,  using  small  smoothing  con¬ 
stants,  produced  the  lowest  error  rate. 
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Conclusion 

The  literature  search  has  revealed  that  the  current 
forecasting  system  used  by  the  SBSS  may  not  provide  the 
best  estimate  of  future  demand.  Additionally,  the  literature 
search  has  shown  the  SBSS  forecasting  method  exhibits  a  sig¬ 
nificant  error  variance  when  compared  to  actual  demand. 

These  facts  lead  to  several  research  questions  which  this 
research  effort  will  explore. 

Research  Questions 

1.  Is  there  a  significant  error  variance  in  the 
multiple  model  forecasting  technique  when  compared  to  pre¬ 
vious  research? 

2.  Which  of  the  forecasting  techniques  selected 
for  this  study  exhibits  the  smallest  error  variance,  i.e., 
which  technique  most  accurately  forecasts  demand? 

Hypothesis 

The  multiple  model  forecasting  technique  performs 
as  well  as,  or  better  than,  single  model  techniques  to  which 
it  is  compared. 
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Chapter  3 


METHODOLOGY 


Introduction 

The  methodology  used  In  this  research  effort 
requires  data  to  be  gathered  and  then  analyzed  using  multi - 
pie  model  forecasting  techniques.  The  results  obtained  from 
the  multiple  model  forecasting  technique  will  be  compared 
to  the  results  of  techniques  used  in  previous  research. 

Using  this  procedure  recommendations  on  the  performance  of 
the  multiple  model  forecasting  technique  can  be  made. 

Model  Selection 

The  current  Standard  Base  Supply  System  (SBSS)  is  a 
mixture  of  moving  averages  and  exponential  smoothing.  It 
is  a  moving  average  since  the  Daily  Demand  Rate  (DDR)  is  a 
quotient  of  cumulative  demands  over  time  divided  by  the 
length  of  time  for  which  these  demands  are  accumulated  (see 
equation  1.7).  It  includes  exponential  smoothing  by  the 
way  updates  to  Cumulative  Recurring  Demand  (CRD)  and  Date 
of  First  Demand  (DOFD)  are  performed  (I9i6).  It  is  not 
apparent,  however,  that  the  current  SBSS  approach  provides 
the  best  estimate  of  future  demand.  An  approach  to  fore¬ 
casting  which  reduces  forecast  error  should  provide  a  more 
accurate  estimate  of  future  demand.  The  techniques 
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considered  in  this  thesis  are  moving  average,  single  expo¬ 
nential  smoothing,  double  exponential  smoothing,  adaptive 
exponential  smoothing,  and  multiple  model  forecasting. 

Moving  Average 

When  demand  for  an  item  does  not  have  a  rapid  growth 
or  seasonal  characteristics,  a  moving  average  can  be  useful 
for  forecasting.  The  moving  average  model  would  be  expected 
to  perform  well  in  normal  constant  demand  situations.  This 
method  assumes  that  the  data  generating  process  constitutes 
a  time  series i 


Fc  "  D  +  e 


(3.1) 


where,  Fc  *  forecasted  demand  for  period  t 

D  =  average  demand  over  time 
e  «  random  error  variable  with  a  mean  of  zero 
as  constant  variance  over  time  (20 i 108) 

The  moving  average  technique  may  be  described  asi 


FC+1  "  °t  1  Dt-1  1  Pt-2  *  *  Dt-n*l  (3.2) 

N 


t-rrH 

1  Y  1 

w  i-t 


where , 


t-rr+1 

I 


i-t 


N 


forecast  for  the  next  period,  t+1 

actual  demand  at  time  t,  t-1,  t-2  ...  t-ntl 

the  number  of  observations  used  in  the 
coverage 
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The  effect  of  moving  average  forecasting  depends  on 
the  number  of  observations  used  in  the  coverage  (N) .  If  N 
is  large,  the  equal  weighting  given  to  each  period  is  small, 
and  random  fluctuations  have  little  effect  on  the  forecast. 

If  N  is  small,  the  model  is  more  sensitive  to  changes  in 
demand  (I2il7). 

Single  Exponential  Smoothing 

In  using  exponential  smoothing,  three  pieces  of  data 
are  needed i  the  most  recent  forecast,  the  actual  demand 
that  occurred  for  that  forecasted  period,  and  a  smoothing 
constant  (Of).  The  equation  for  a  single  exponential  smooth¬ 
ing  forecast  isi 

f»l  -  Ft  ♦  ar(Dt  -  Ft)  (3.3) 

where,  F.. ,  ■  the  exponentially  smoothed  forecast  for 
m  period  Vl 

Ft  ■  the  exponentially  smoothed  forecast  for  the 
prior  period  t 

Dt  ■  the  actual  demand  in  the  prior  period  t 

CK  ■  the  desired  response  rate,  or  smoothing 

constant  (I3i58) 

This  equation  states  that  the  new  forecast  is  equal  to  the 
old  forecast  plus  an  adjustment  proportional  to  the  differ¬ 
ence  between  the  previous  forecast  and  the  actual  experience. 
The  closer  Ct  is  to  1,  the  more  the  new  forecast  will  incor¬ 
porate  an  adjustment  for  the  error  in  the  upcoming  forecast. 
The  closer  Qf  is  to  0,  the  less  sensitive  the  new  forecast 
will  be  to  the  error  in  the  prior  forecast  (I3i58). 
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Double  Exponential  Smoothing 


Double  exponential  smoothing  is  an  extension  of 
single  exponential  smoothing.  It  is  usually  applied  to 
items  that  exhibit  a  trend  pattern.  Both  single  and  double 
smoothed  values  lag  actual  data  when  a  trend  exists,  the 
difference  between  single  and  double  smoothed  values  can  be 
added  to  the  single  smoothed  values  and  thereby  adjust  for 
trend  (I9ill).  The  single  exponential  model,  equation  (3.3), 
can  be  rewritten  asi 


St  “  St-1  +  « (D t  “  St-l)  (3.4) 

where,  "  exponentially  smoothed  demand  in  period  i. 

The  apparent  trend  component  of  the  error  factor  of  our 
forecast  isi 


« 


-  S 


t-1 


(3.5) 


This  trend  component  can  also  be  smoothed  or  averaged  over 
time  usings 

Tt  *  Tt-i  ♦  a(ij  -  it.t)  (3.6) 

The  double  exponentially  smoothed  forecast  or  the  smoothed 
average  forecast  including  the  trend  component  is  the  simple 
smoothed  average  plus  the  smoothed  trend  component  as  cor¬ 
rected  for  lag  in  Tt  by  the  term  (1 


t+1 


S  +  — =— T 

t  a  t 


(3.7) 
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In  order  to  use  this  approach  only  three  data  values 
and  a  smoothing  constant  are  required.  Use  of  this  model 
should  be  based  on  assumptions  about  the  demand  pattern. 

If  the  time  series  has  a  trending  average  demand  rate,  the 
double  exponential  smoothing  model  produces  an  accurate 
estimate  of  demand  (I2il9). 


Adaptive  Exponential  Smoothing 

The  adaptive  smoothing  techniques  adjust  Qf  over 
time  based  on  the  size  of  the  forecasting  error.  The  adap¬ 
tive  model  utilized  here  is  the  same  as  in  single  exponen¬ 
tial  smoothing  with  one  exception.  The  value  for  (X  is 
derived  from  the  equationi 


a 


t+i 


M. 


(3.8) 


where,  Et  ■  /?(et)  + 
Mc  -  /3|«t|  ♦ 


D  -  F 
t  t  t 

P  ■  smoothing  constant 
and  |  |  denotes  absolute  values  (16i54). 

The  characteristic  of  not  specifying  a  value  for  CX 
is  attractive  when  thousands  of  items  require  forecasting. 
Also,  this  method  can  change  the  value  of  Of  when  changes 
in  the  pattern  of  the  data  have  made  the  initial  Of  value 
no  longer  appropriate  (16i53). 
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Multiple  Model  Forecasting 

The  model  that  most  accurately  forecasts  demand  will 
be  selected  from  all  models  considered.  In  order  to  do 
this,  each  demand  history  will  be  broken  into  three  periods i 
a  base  (or  start -up)  period,  a  test  period,  and  a  prediction 
period. 

The  base  (or  start-up)  period  is  the  historical  data 
base  required  for  each  model.  The  length  of  the  base  period 
varies  from  model  to  model. 

After  the  forecasting  start  point  is  established, 
each  model  forecasts  demand  for  a  given  test  period.  The 
forecasts  are  compared  to  actual  demand  for  the  test  period. 
The  model  with  the  smallest  variation  from  the  actual  demand 
is  selected  to  forecast  demand  for  the  next  prediction 
period. 

This  is  a  recurring  process  in  that  the  forecasting 
horizon  defines  a  new  prediction  periodi  the  past  prediction 
period  becomes  the  current  test  period i  the  past  test  period 
joins  the  base  period  data  bases  and  older  data  are  dropped 
from  the  base  period  (13il2). 

The  multiple  model  method  employed  in  this  research 
is  shown  in  Appendix  D,  subroutine  "Multi ."  The  actual 
demand  for  the  four  previous  forecasting  periods  is  added 
and  stored  as  a  new  variable.  The  values  forecasted  for 
these  same  four  periods  by  each  technique  are  added  and  this 
total  is  subtracted  from  the  total  for  the  actual  demand. 
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This  process  is  repeated  for  each  forecasting  technique. 

The  technique  which  displays  the  lowest  variance  from  the 
actual  demand  is  selected  by  the  multiple  model  technique 
as  the  technique  it  will  use  to  forecast  demand  for  the 
next  period. 

Accuracy  of  Forecasting  Models 

Demand  for  an  item  is  generated  through  the  inter¬ 
action  of  a  number  of  factors.  This  interaction  is  extremely 
complex.  Due  to  this  complexity,  all  forecasts  will  contain 
some  error. 

In  many  forecasting  situations  accuracy  is  treated 
as  the  main  criterion  for  selecting  a  forecasting  method. 

In  spite  of  this  fact,  little  systematic  work  has  been  done 
to  develop  a  framework  for  measuring  and  evaluating  accuracy- 
related  issues.  One  of  the  difficulties  in  dealing  with  the 
criterion  of  accuracy  in  forecasting  situations  is  the 
absence  of  a  single  universally  accepted  measure  of  accuracy 
(16i569).  When  discussing  forecast  errors  it  is  appropriate 
to  distinguish  between  sources  of  error  and  the  measurement 
of  error. 

Sources  of  error  can  be  classified  as  either  bias 
or  random.  Bias  errors  occur  when  a  consistent  "mistake” 
is  made.  Random  errors  can  be  defined  as  those  that  cannot 
be  explained  by  the  forecast  model  employed. 
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Several  of  the  common  terms  used  to  measure  the 
degree  of  error  are  mean  absolute  percentage  error  and  mean 
absolute  deviation.  Additionally,  measures  of  bias  may  be 
used  to  indicate  the  amount  of  any  positive  or  negative  bias 
in  the  forecast  (9i24l), 


Mean  Absolute  Percentage  Error 

One  approach  used  to  detect  forecast  error  is  to 
take  the  absolute  value  of  each  individual  percentage  error 
to  obtain  the  mean  absolute  percentage  error  (MAPE) . 

Percentage  error  is  a  relative  measure  of  accuracy 
and  is  defined  asi 


PEt  - 


D  -  F. 


(100) 


(3.9) 


where,  PEt  ■  percentage  error  for  any  time  period  t 
Dt  -  actual  value  for  time  period  t 
Fc  ■  forecast  value  for  time  period  t  (I6i570) 

The  equation  for  computing  the  MAPE  over  N  time 


periods  is  given  byi 


MAPE 


N 


(3.10) 


where,  MAPE  *  mean  absolute  percentage  error 

PE^  ■  percentage  error  for  a  time  period  (I5i570) 

However,  the  MAPE  Itself  does  not  always  give  a  good 

basis  of  comparison  as  to  the  gains  in  accuracy  made  by 
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applying  a  specific  forecasting  method.  MAPE  indicates  how 
well  the  forecasting  model  fits  the  data.  In  this  way,  it 
aids  in  evaluating  how  accurate  the  model  is  (16i568). 

Mean  Absolute  Deviation 

The  most  straightforward  measure  of  error  is  mean 
absolute  deviation  (MAD) .  It  is  computed  using  the  differ¬ 
ences  between  the  actual  demand  and  the  forecasted  demand 
without  regard  to  sign.  Since  it  is  the  mean  deviation,  the 
sum  of  the  absolute  deviations  is  divided  by  the  number  of 
data  points.  Stated  in  equation  formt 

I  |Dt  - 

MAD  »  -  (3.11) 

N 

where,  »  actual  demand  for  period  t 

Ft  ■  forecasted  demand  for  period  t 
N  ■  total  number  of  periods 

and  |  |  denotes  the  absolute  value  (9i242). 

The  mean  absolute  deviation  ignores  whether  the  forecast  is 

greater  or  less  than  actual  demand  and  simply  measures  the 

magnitude  of  difference. 

Bias 

Bias  indicates  the  directional  tendency  of  forecast 
errors.  If  the  forecasting  technique  used  consistently  over¬ 
estimates  actual  demand,  bias  will  be  positive*  consistently 
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cause  the  bias  to  be 


underestimating  actual  demands  will 
negative.  Bias  is  defined  as< 

n 

£  (Ft  -  Vi 

a  -j”! _ _ 


where,  Fc  »  forecasted  demand  for  period  t 
Dc  “  actual  demand  for  period  t 
N  »  total  number  of  periods  (Ii33i) 

Conclusion 

All  three  techniques  described  here  are  used  to 
evaluate  multiple  model  forecasting.  The  measures  of  error 
in  the  forecasting  techniques  are  used  to  determine  which 
forecasting  method  will  provide  the  most  accurate  estimate 
of  demand.  By  using  the  forecasting  technique  which  is  the 
most  accurate,  the  error  variance  in  the  Standard  Base  Supply 
System  can  hopefully  be  reduced. 
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Chapter  4 


DATA  COLLECTION  AND  ANALYSIS 


Introduction 

The  purpose  of  this  chapter  Is  to  Identify  the 
source  of  data  used  In  the  analysis  and  describe  the  selec¬ 
tion  procedure  used  to  obtain  a  sample  of  279  National 
Stock  Numbers  (NSNs).  Additionally,  this  chapter  will  pre¬ 
sent  an  analysis  of  that  data. 

Data  Collection 

As  mentioned  in  Chapter  1,  this  research  had  the 
sponsorship  of  the  Air  Force  Logistics  Management  Center 
(AFLMC),  Gunter  AFS,  Alabama.  Two  computer  tapes  were  pro¬ 
vided  which  contained  Standard  Base  Supply  System  data  from 
Dover  AFB,  Delaware,  for  a  two  and  one-half  year  period 
fro.n  October  1976  to  March  1979.  One  tape  contained  item 
record  data  for  33  stock  classes,  ranging  from  the  1005  to 
the  5315  stock  class.  The  total  number  of  stock  classes  was 
115  containing  a  total  of  10,607  individual  National  Stock 
Numbers.  The  remaining  tape  contained  transaction  history 
data  for  these  same  NSNs.  A  Fortran  program  (Appendix  A) 
was  developed  which  aggregated  these  records  by  stock  class) 
total  items  in  each  class t  and  total  recurring  demands. 
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The  next  step  was  selecting  the  stock  classes  to  be 
used  in  the  actual  analysis.  A  previous  study  by  Patterson 
(19)  testing  alternative  forecasting  techniques  for  the 
Standard  Base  Supply  System  (SBSS)  used  stock  classes  28,  47, 
59,  and  66.  It  was  the  researchers'  intention  to  select  a 
sampling  of  these  same  classes  in  order  to  make  a  comparative 
analysis  of  results.  However,  since  the  data  only  Included 
through  the  53  class,  two  additional  alternative  stock 
classes  were  chosen.  A  Fortran  program  (Appendix  B)  was 
used  to  prioritize  the  data  into  groups  of  high  activity 
(frequency  of  demand)  and  quantity  demanded.  With  this 
information,  the  authors  were  able  to  identify  the  51  and 
53  stock  classes  as  having  characteristics  that  would  lend 
themselves  to  alternative  forecasting  methods  and,  partic¬ 
ularly,  to  the  multiple  model  method.  Table  1  shows  a  cap- 
sulized  summary  of  the  activity  and  demands  for  the  four 
stock  classes  selected. 


Table  1 

Profile  of  Preliminary  Data 


Stock  Class 

Total  NSN's 

Total 

Recurring  Demands 

28 

386 

18,196 

47 

1313 

84,939 

51 

2407 

126,571 

53 

2331 

842,654 
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Due  to  the  large  frequency  and  quantity  of  demand 
data  and  the  limited  computer  resources,  a  sampling  of  this 
population  was  taken.  Recurring  demands  were  tabulated  from 
the  transaction  history  tape  (see  Fortran  program  in 
Appendix  C)  and  aggregated  into  130  weeks  or  data  points  for 
each  selected  item. 

Prior  to  the  selection  of  the  final  sample,  the  data 
were  scanned  for  any  obvious  outliers.  As  in  Patterson's 
research,  one  stock  number  (47  stock  class)  had  6,160  units 
requested  in  one  period  and  zero  demands  for  the  remaining 
129  periods.  This  item  was  dropped  from  the  sample.  A  28 
stock  class  item  was  also  discarded  because  data  were  found 
to  be  available  only  for  the  first  one  and  one-half  years  of 
the  survey  period. 

Table  2  is  a  summary  of  the  final  stock  classes 
used  in  testing  the  alternative  forecasting  methods. 

Table  2 

Profile  of  Stock  Classes  Selected  for  Analysis 


Stock  Class 

Total  NSN's 

Total 

Recurring  Demands 

28 

99 

8,137 

47 

99 

14,131 

51 

31 

21,659 

53 

50 

234,689 
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Summary 


In  summary,  the  data  used  for  the  actual  analysis 
consisted  of  centrally  procured  items  identified  by  the 
Expendability,  Repairability,  Recoverability  Category  (ERRC) 
of  XB3.  A  total  sample  of  279  line  items  was  selected 
based  on  high  activity  and  a  relatively  high  frequency  of 
reorder  periods.  An  inherent  bias  in  the  data  collection 
was  the  researchers*  attempt  to  identify  those  stock  numbers 
which  exhibited  characteristics  which  would  accommodate  the 
time  series  forecasting  techniques  to  be  exercised. 

This  research  addressed  the  applicability  of  multi¬ 
ple  model  forecasting  as  an  alternative  to  the  standard  base 
supply  forecasting  system.  In  order  to  test  this  hypothesis, 
the  data  collected  by  the  methods  described  were  processed 
by  use  of  a  Fortran  IV  computer  program  (see  Appendix  D) . 

Computer  Output 

For  each  stock  number  in  each  stock  class,  the 
actual  units  demanded,  the  units  forecasted  by  each  fore¬ 
casting  technique  used,  and  the  units  forecasted  by  the  mul¬ 
tiple  model  forecasting  technique  were  shown  on  the  output 
(see  Appendix  E) .  The  mean  absolute  deviation  (MAO)  com¬ 
paring  the  forecasted  demand  and  the  actual  demand  was  com¬ 
puted  (see  equation  3.11).  Additionally,  a  bias  (equation 
3.12)  and  mean  absolute  percentage  error  (MAPE,  equation 
3.10)  were  computed.  These  values  were  necessary  to  complete 
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the  statistical  tests  required  to  determine  the  accuracy  of 
the  various  forecasting  methods. 


Data  Base  Comparison 

In  order  to  attain  the  objective  of  evaluating  the 
multiple  model  forecasting  method,  it  is  necessary  to  show 
that  a  similar  data  base  exists  with  which  to  compare  our 
results.  In  his  research,  Patterson  obtained  a  mean  abso¬ 
lute  deviation  for  single  and  double  exponential  smoothing, 
using  an  alpha  of  .2,  among  others  (19i9).  A  comparison  of 
the  means  for  these  two  forecasting  techniques  will  test  the 
data  base  used  to  see  if  it  is  essentially  the  same  as  that 
used  by  Patterson  in  his  research.  The  precise  stock  numbers 
used  by  Patterson  were  unknown  to  us. 

To  test  this  assumption,  95  percent  confidence  inter¬ 
vals  were  computed.  It  is  necessary  to  calculate  a  value  for 
the  standard  deviation  for  Patterson's  research,  since  none 
was  given.  When  forecast  errors  are  normally  distributed, 
the  MAD  is  related  to  the  standard  deviation  byi 


1  Standard  Deviation 


where,  7 T  "  3.1416 
or  conversely. 


(4.1) 


1  MAD  ■  0.8  Standard  Deviation  (9i242)  (4.2) 
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The  central  limit  theorem  tells  us  that  for  almost  all  popu¬ 
lations  the  sampling  distribution  of  the  sample  mean  is 
approximately  normal  when  the  sample  size  is  sufficiently 
large  (I3i202).  Since  the  sample  size  used  by  Patterson  is 
large  (200),  normality  can  be  assumed,  and  a  value  can  be 
calculated  for  the  standard  deviation. 

Using  this  value  for  the  standard  deviation  of  both 
samples,  a  95  percent  confidence  interval  was  constructed. 

To  construct  a  confidence  interval,  an  unbiased  estimator  of 
the  difference  in  the  two  sample  means  is  required.  This 
estimator  is  denoted  by  Di 

D  -  7  -  X  (4.3) 


where,  Y  ■  mean  of  observations  from  sample  2 

X  *  mean  of  observations  from  sample  1  (I8i3l2) 

Additionally,  an  unbiased  estimator  for  the  standard  devia¬ 
tion  of  D  is  required i 


S2(D)  ■  S2(Y)  ♦  S2(X) 


(4.4) 


where,  S2(Y)  *  S22 


n 


S2(X)  -  (18 1 3.6) 


The  formula  for  calculating  the  two-sided  confidence  inter¬ 
val  1st 


L  -  5  -  3(1-CK/2)S(D) 

U  -  B  +  3(1-072 )S(U)  (18.316) 


(4.5) 
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The  results  of  the  confidence  Interval  calculations 


are  shown  In  Table  3: 


Table  3 

95%  Confidence  Interval  for  Data  Base  Comparison 


Model 

Lower 

Limit 

Mean 

MAD 

Upper 

Limit 

Stock  Class  23 

Single  Exponential 

.6282 

2.11 

2.5918 

Double  Exponential 

.8670 

2.36 

2.8931 

Stock  Class  47 

Single  Exponential 

.5062 

2.06 

2.2938 

Double  Exponential 

.2405 

2.47 

2.6396 

Since  the  means  fell  within  the  limits*  we  have  concluded 
that  for  practical  purposes  our  data  base  and  that  used  by 
Patterson  were  essentially  the  same. 

Data  Analysis 

The  results  of  the  computer  analysis  of  the  data 
are  shown  in  Table  4,  In  three  of  the  four  stock  classes* 
on  the  average*  single  exponential  smoothing  consistently 
displayed  the  lowest  mean  absolute  deviation.  This  is  con¬ 
trary  to  what  was  expected.  The  data  were  selected  from 
the  various  stock  classes  based  on  large  amounts  of  demand 
occurring  in  the  most  time  periods.  Essentially  the  "deck 
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Table  4 


Results  of  Computer  Analysis 


Mean  MAD 

Mean  Bias 

Mean  MAPE 

Stock  Class  28 

Single  Exponential 

2.1115 

• 

.6691 

88.7535 

Double  Exponential 

2.3659 

- 

.3165 

114.7806 

Adaptive  Exponential 

2,9618 

.5948 

178.4009 

4 -Term 

2.3750 

- 

.2396 

117.1407 

8 -Term 

2.2967 

- 

.3491 

111.1047 

Multiple  Model 

2.7200 

• 

.2081 

152.5846 

Stock  Class  47 

Single  Exponential 

2.0682 

• 

1.0852 

213.0647 

Double  Exponential 

2.4756 

- 

.1933 

278.1608 

Adaptive  Exponential 

3.8648 

1.3457 

388.1105 

4-Terra 

2.5769 

• 

.2905 

265.9312 

8-Term 

2.3863 

- 

.4338 

251.1996 

Multiple  Model 

2.8990 

.6206 

328.0684 

Stock  Class  51 

Single  Exponential 

13.6679 

m 

1.6520 

615.5107 

Double  Exponential 

14.8354 

.8196 

781.0071 

Adaptive  Exponential 

18.1501 

3.8259 

1081.6932 

4 -Terra 

14.4877 

.3390 

710.0943 

8-Term 

13.9992 

.6066 

677.3004 

Multiple  Model 

16.1558 

'1.8196 

894.6500 

Stock  Class  53 

Single  Exponential 

98.2697 

m 

3,9447 

4023.2286 

Double  Exponential 

107.5422 

14.2225 

5036.5975 

Adaptive  Exponential 

114.3334 

15.5944 

5373.4977 

4 -Term 

103.0422 

m 

1.0954 

4245.1873 

8 -Terra 

97.7790 

m 

2.0539 

4073.0743 

Multiple  Model 

107.2436 

8.0751 

4762.1244 

36 
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was  stacked"  In  favor  of  the  time  series  techniques.  The 
multiple  model  technique  did  not  produce  the  most  accurate 
estimate  of  demand,  nor  did  it  even  tie  with  the  most  accu¬ 
rate  method.  In  the  case  of  the  51  and  53  stock  classes, 
the  MAD  for  all  techniques  considered  rises  considerably. 

This  indicates  the  problems  encountered  in  using  time  series 
analysis  on  data  that  are  highly  erratic  and  sparse  (many 
zero  entires).  The  MAPE  reflects  this  highly  erratic 
nature.  The  relatively  large  values  for  the  MAPEs  indicate 
abrupt  changes  in  the  demands  in  the  data  ( see  Appendix  E) . 
There  are  numerous  periods  in  which  zero  demand  is  incurred 
and  suddenly  a  large  change  in  demand  occurs  and  abruptly 
drops  off  to  zero  again. 

Conclusion 

Through  the  use  of  confidence  intervals  it  was  shown 
that  the  data  base  used  in  this  research  was  essentially  the 
sane  as  that  used  by  Patterson  in  his  previous  research.  The 
Fortran  IV  program  produced  results  that  indicated  single 
exponential,  and  not  multiple  model  forecasting,  was  the 
most  accurate  of  the  techniques  considered.  The  implica¬ 
tions  and  conclusions  which  have  been  drawn  from  this  will 
be  further  discussed  in  the  next  chapter. 
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Chapter  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


Introduction 

This  research  was  directed  at  exploring  the  hypoth¬ 
esis  i  the  multiple  model  forecasting  technique  performs  as 
well  as,  or  better  than,  single  model  techniques  to  which 
it  is  compared.  The  research  addressing  this  hypothesis  led 
to  several  conclusions,  which  can  be  associated  with  the 
four  objectives  of  this  thesis. 

Conclusions 

The  first  objective  of  this  research  was  to  evaluate 
the  multiple  model  forecasting  technique  as  an  alternative 
to  the  Standard  Base  Supply  System  forecasting  technique. 
Closely  related  to  this  was  the  second  objective,  to  deter¬ 
mine  which  forecasting  approach  will  provide  the  most  accu¬ 
rate  estimate  of  future  demand. 

The  multiple  model  forecasting  technique  and  the 
additional  techniques  were  processed  by  the  use  of  the  For¬ 
tran  IV  program  referred  to  in  Chapter  4.  Table  4  in  Chap¬ 
ter  4  reflects  the  results  of  the  data  analysis.  As  can  be 
seen,  multiple  model  forecasting  did  not  produce  the  most 
accurate  estimate  of  demand.  Based  on  the  mean  MAD  computed 
after  the  computer  analysis  was  completed,  multiple  model 
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consistently  ranked  next  to  last  in  accuracy.  On  three  out 
of  the  four  stock  classes,  single  exponential  smoothing  was 
ranked  as  the  most  accurate.  On  the  fourth  class,  stock 
class  53,  single  exponential  smoothing  was  the  second  most 
accurate  method,  following  closely  behind  8 -terra  moving 
average. 

The  bias  is  an  indicator  of  the  directional  tendency 
of  the  forecast  errors  for  each  forecasting  technique.  In  a 
perfect  forecasting  model  the  bias  would  be  expected  to  be 
zero,  i.e.,  the  random  errors  due  to  overforecasting  would 
be  offset  by  random  errors  due  to  underforecasting.  Again, 
referring  to  Table  4,  in  three  out  of  the  four  stock  classes 
analyzed,  multiple  model  forecasting  consistently  on  the 
average  overforecasts  the  demand.  The  Implications  of 
overforecasting  demand  come  into  play  when  the  holding  cost 
of  an  item  is  considered.  Holding  cost  is  important  to  any 
manager  concerned  with  any  form  of  inventory.  Holding  cost, 
along  with  order  costs  and  demand  rates,  makes  up  the  eco¬ 
nomic  order  quantity  (EOQ) .  The  multiple  model  forecasting 
technique  used  in  this  research  would  cause  the  holding  cost 
to  rise,  if  only  slightly,  due  to  its  consistent  over¬ 
estimating  of  demand  and  subsequent  over- stocking  of  inven¬ 
tory.  Single  exponential  smoothing,  as  used  in  this 
research,  consistently  tinder-estimates  demand  and  therefore 
again  becomes  a  more  acceptable  forecasting  technique  when 
compared  to  others,  including  multiple  model. 


Finally,  the  mean  absolute  percentage  error  (MAPE) 
of  the  multiple  model  technique  consistently  ranks  as  the 
next  to  least  accurate.  The  MAPE  indicates  how  well  the 
technique  in  question  fits  the  data.  Therefore,  the  multi¬ 
ple  model  forecasting  used  in  this  research  does  not  fit 
the  data  as  well  as  other  methods.  Single  exponential 
smoothing  fits  the  data  best,  i.e.,  it  consistently  exhib¬ 
ited  the  lowest  MAPE.  The  high  MAPEs  encountered  in  this 
research  indicate  the  highly  erratic  nature  of  the  SBSS. 

An  example  of  the  erratic  nature  of  the  demand  data  is  shown 
in  Appendix  E.  Additionally,  It  can  be  seen  that,  as  a  gen¬ 
eral  rule,  trend  is  not  apparent  in  SBSS  items.  Therefore, 
methods  of  forecasting  that  did  adjust  for  trend  did  not 
reduce  the  forecast  error.  When  an  abrupt  change  occurred 
in  the  data,  both  double  and  adaptive  exponential  smoothing 
reacted  to  the  change  as  if  it  were  a  trend.  However,  the 
data  usually  abruptly  returned  to  zero  and  a  trend  did  not 
develop.  As  a  result,  these  forecasting  methods  were  not 
among  the  most  accurate  of  the  methods  tested. 

The  third  objective  of  this  research  was  to  deter¬ 
mine  the  effectiveness  of  implementing  the  proposed  system. 
In  order  to  demonstrate  a  significant  cost  savings  to  the 
Air  Force,  there  must  be  a  significant  increase  in  fore¬ 
casting  accuracy.  The  results  indicate  that  the  multiple 
model  forecasting  technique  used  here  would  not  support 
this  significant  increase  in  forecasting  accuracy. 


Therefore,  any  cost  savings  attributable  to  implementing  the 
multiple  model  forecasting  technique  are  in  doubt.  It  also 
becomes  apparent  that  the  time  series  forecasting  models 
used  in  this  research  did  not  yield  any  significant  improve¬ 
ment  in  the  forecasting  of  demand  for  SBSS  items. 

The  final  objective  of  this  research  was  to  recom¬ 
mend  actions  based  upon  the  results  obtained.  The  following 
section  will  address  this  objective. 

Recommendati ons 

Implied  throughout  this  study  is  the  fact  that 
inventory  models  are  used  to  assist  managers  in  determining 
when  items  should  be  ordered  and  how  large  the  order  should 
be  so  as  to  minimize  variable  costs  and  insure  that  adequate 
stockage  is  on-hand  to  meet  mission  requirements.  Using  the 
EOQ  model  as  a  depth  model,  however,  has  two  inherent 
assumptions,  i.e.,  the  demand  activity  remains  relatively 
continuous  and  the  items  are  relatively  low  in  cost.  Refer¬ 
ring  to  Table  5,  we  note  that  this  model  was  intended  for 
use  under  the  conditions  which  are  present  in  Quadrant  I. 

Table  5 

Inventory  Characteristics 


Inventory  Item  Cost 
HI  LO 


HI 

II 

I 

I  tern 

Demand 

LO 

III 

IV 
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Quadrant  I  exhibits  the  ideal  conditions  for  EOQ,  l.e,,  a 
high  frequency  of  demand  with  continuous  demand  behavior. 
Scarcity  of  resources  and  other  fiscal  environmental  factors 
have  made  the  conditions  within  Quadrant  II  not  at  all  un¬ 
common.  Average  item  costs  have  been  increasing  signifi¬ 
cantly  in  recent  years.  The  EOQ  model,  while  not  optimal, 
can  handle  these  conditions  fairly  well  as  long  as  the 
demand  patterns  remain  relatively  continuous.  In  addition 
to  costs,  technology  factors  have  increased  reliability  on 
many  items  resulting  in  lower  demand  frequency.  This  has 
the  effect  of  moving  many  items  into  Quadrant  III  and  IV. 

An  example  would  be  vacuum  tubes  which  are  now  replaced  by 
solid  state  technology  resulting  in  a  higher  Mean  Time  Be¬ 
tween  Failure.  Quadrant  III  and  IV  contain  characteristics 
of  low  demand  for  both  low  and  high  value  items.  The  math¬ 
ematical  assumptions  inherent  within  the  EOQ  model  were  not 
meant  to  address  the  conditions  of  low  demand,  especially 
the  sparse  lumpy  conditions  found  frequently  in  the  SBSS 
data.  Examples  of  this  demand  data  are  found  in  Appendix  E. 

We  conclude  that  the  EOQ  model  is  an  inappropriate 
depth  model  for  items  characterized  by  erratic  and  sparse 
demand  patterns.  Our  recommendation  is  to  divert  further 
efforts  from  Improving  SBSS  forecasting  techniques  to  direct¬ 
ing  research  efforts  in  a  search  for  alternative  decision 
models  which  will  be  more  appropriate  for  the  SBSS  depth  of 
stockage  decision. 
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APPENDICES 
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APPENDIX  A 

DATA  RETRIEVAL  PROGRAM 
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10U  CHARACTER  STCLASS*4,CLHOLD*4.STNUMB*15,SNHOLD*1 5 

20*  INTEGER  CLASSCNT.TOTCNT.RDCNT.RD 
•JO#  READC20, 100,END=900)STCLASS,STNUMB,RD 

70#1 00  F0RNAT<T5,A4.T5.A15.T116,I6) 

80#  T0TCNT»1 

90#  CLASSCNT=1 
100#  CLHOLOaSTCLASS 

110#  SNHOLD-STNUMB 

120#  RDCNT=RD 

130#  ITEMCNT=1 

140#  URITE(6, 125) 

15014125  FORNATCMSTOCK  CLASS  TOTAL  ITEMS  TOTAL  RECURRING  DEHANB") 

160#  DO  500  K«1, 11000 

170#  READC20, 100,  END=900)STf;_  ASS,  S  THUMB,  RD 

180#  T0TCNT=TOTC#T+1 

220#  IF ( STNUMB . GE . SNHOLD ’ GO  TO  150 

230#  PRINT,"  ST  NUMBERS  OUT  OF  SEQsREC#  ".TOTCNT 

240#  STOP 

250#1 50  SNHOLD=STNUMB 

260#  IF(STCLASS.EQ.CLHOLD)GO  TO  200 

330#  CLASSCNT»CLASSCNT+1 

340#  URITE<6, 1010)CLH0LD,ITEMCNT,RDCNT 

350#  CLHOLD-STCLASS 

360#  RDCNT=0 

370#  ITEMCNT*0 

380H200  RDC#T*RDCNT4RD 

390#  1TEMCNT=ITEMCNT+1 

400#500  CONTINUE 

410#900  URITE16, 1010JCLH0LD, ITEMCNT,RDCNT 
420#  URITE(6,950)CLASSCNT 

430#950  FORMATC///-  TOTAL  NUMBER  OF  STOCK  CLASSES  *  -,I6) 

440#  URI TEC  6, 1000) TOTCNT 

450R1000  FORMAT!"  TOTAL  NUMBER  OF  RECORDS  »  ”,I4) 

460  # 1010  F0RHAT(/2X,A4,T18,I?fT33,I11> 

470#  STOP 

480#  END 
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APPENDIX  B 

SELECTION  PROGRAM  FOR  RANKING  BY 
DEMAND  AND  FREQUENCY 
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LIST  GE09 


10  INTEGER  SC,7D, IARRAY (131) 

20  CHARACTER*^  SN.HOLDSN 
30  IFLA6=0 

40  HOL DSN3’ 123454789 123456“ 

50  50  READ140, 100. END=500)SC,SN. IDATE. IDEMAND 
60  100  FORMAT <T1 ,12,T1,A15,T67,I4,T79,I6) 

70  IF(SC.6E.S3)G0  TO  200 
VO  GO  TO  50 

100  200  IF (SN.NE.HOLDSN)GO  TO  300 

110  CALL  TALLY ( SN, IDATE , IARRAY , I DEMAND , ITOT ) 

120  GO  TO  50 

130  300  IF ( IFLAG.EQ.O)GO  TO  350 

140  DO  400  K=1,131 

150  IBUCKETS=IBUCK£TS+IARRAY(K) 

160  400  CONTINUE 

170  IBUCKETSsIBUCKETS+ 100000 

180  I TOT3ITQT+ 1000000000 

190  URITE(50,425)H0LDSN, IBUCKETS,ITOT 

200  425  FORMAT (1X,At5,1X,14,1X,I8) 

205  350  IF(SC.GT.53)G0  TO  500 

210  IFLAG31 

220  HOLDSN=SN 

225  IBUCKETS=0 

230  IT0T30 

240  DO  450  J=1 ,131 

250  IARRAY1 J)=0 

260  450  CONTINUE 

270  CALL  T ALL Y ( SN, IDATE, I ARRAY, IDEMAND, ITOT) 

280  GO  TO  50 

320  500  URITE<6,550) 

330  550  FORMAT <"  END  OF  SN  PROGRAM") 

340  STOP 
350  END 

360  SUBROUTINE  TALLY1SN. IDATE, IARRAY, IDEMAND, ITOT) 
370  CHARACTER  SN*15 
380  INTEGER  IARRAY<131) 

390  IF1IDATE. GE. 4276. AND. IDATE. LE. 9091 )G0  TO  600 
400  URITE(6,575)SN, IDATE 
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410  57S  F0RHAT1"  DATE  OUT  OF  BOUNDS*  *.At5,5X>I4) 

420  RETURN 
430  600  CONTINUE 

440  IF < IDATE. G£. 6276. AND. IDA7E.LE. 6366) I SUB=( I DATE-626?)/? 

450  IF ( IDATE.GE. 7001 .AND. I BATE. LE. 7365) ISUB=( I DATE-6994 )/?+1 3 

460  IF ( IDATE.GE .8001 .AND. IDATE. LE.8365)ISUB=(IDATE-7993)/7+65 

470  IF (IDATE.GE. 9001. AND. IDATE. LE. 9091 )ISUB=(IDATE-8992>/7+11? 

480  I ARRAY < ISUB)  =  1 

490  IT0T=ITO7+IBEMANB 

500  RETURN 

510  END 


LIST  S0RT-1B 
0099INN0RN, J 

01004: I DENT :UP1 186, AFIT  ACDS  LT  NAIER 

01204:GNAP:NDECK 

013O:600SH:N0HAP 

0140:S0RT:IN0UT,,5 

0150:FIELD: (C16,C5,C9) 

0160:SEQ:(D2.D3,A1 ) 

0170:FILCB:INQUT.**,2 

0180:END 

01904:EXECUTE 

02004:LINITS:10 

02 1 04 :PRNFL : SA , R , S . BOBF /DAT  A 1 -47 

02204 :PRHFL:SZ,U,S.B0BF/DATA1 -B 

02304: FILE: SI ,X1R,10R 

02404 : F ILE :S2,X2R, 1 0R 

02504 : FILE :S3, X3R , 1 OR 

02604:ENDJ0B 


LIST  SORT-ID 
0099N»NORH,J 

01004:IDENT:UF1 186, AFIT  ACDS  LT  NAIER 

01204:GHAP:NDECK 

0130:600SN: NONAP 

01 40: SORT : INOUT ,  ,5 

01  SO: FIELD: (C16.C5.C9) 

01 60:SE0: (D3.D2.A1 ) 
0170:FILCB:IN0UT,**,2 
0180:END 
01904 :EXECUTE 
0200*  sLIHITS :  1 0 

02104:PRNFL:SA.R,S,B0BF/DATA1 -47 
02204:PRHFL:SZ,U,S.B0BF/DATA1-D 
02304:FILE:S1 .X1R.10R 
02404:FILE :S2,X2R. 1 0R 
02504 : FILE : S3, X3R,1 OR 
02604 :ENDJOB 
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APPENDIX  C 

PROGRAM  FOR  DATA  ANALYSIS 
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LIST  NO VEIT 


10  CHARACTER*^  SN1 ,SN2 

20  INTEGER  BUKET1 < 1 30) ,8UKET2< 130) , BUKET31 130) 

30  5  READ( 10, 1 0,END=500)SN1 , (BUKET1 (K) ,K=1 , 130 ) 

40  READ(20, 1 0,END=600)SN2, ( BUKET2 ( J ) , J=1 , 1 30) 

50  10  FO RNAT<A15,5(/26I5>> 

60  IF ( SN1 .NE.SN2 >G0  TO  700 
70  DO  30  1=1,130 

80  BUKET3(L)*BUKET1 <L)+BUKET2(L> 

90  30  CONTINUE 

100  URITE(30,50)SN1,(BUKET3< I >,1=1,130) 

110  50  FORNAT (A1 5, 12I5,7(/15I5)/13I5) 

120  GO  TO  5 

130  500  URITE(6,550 > 

140  550  FORHAT ("  END  OF  FILE-1") 

150  READ(20,10,END=575) 

160  WRITE < 6 , 590 )SN2 

170  590  FORNAK"  E0F-1,BUT  INFO  STILL  ON  FILE-2. SN2=" ,A15> 
180  STOP 

190  575  URITE(6,595) 

200  595  FORNAT ( "  EOF  ON  BOTH  FILES,  ALL  OK") 

210  STOP 

220  600  URITE(6,660)SN1 

230  660  FORNAT  1 "  UNEXPECTED  E0F-2,SN1  =  '',AI5) 

240  STOP 

244  700  URITE(6,750)SN1 ,SN2 

246  750  FORHAT ( "  UNMATCHED  SN'S=".A15,2X,A15) 

248  STOP 
250  END 

410  I ARRAY ( ISUB  >=I ARRAY (ISUB )  +  IDENAND 

420  IF ( I ARRAY ( ISUB) . GT .9999?)URITE<6,700)SN 

430  700  FORNAT <  "IDENAND  HAS  EXCEEDED  99999,  SN=",A15) 

440  RETURN 
450  END 

460  SUBROUTINE  INIT ( IFOUND, I ARRAY ) 

470  INTEGER  I ARRAY (130) 

480  DO  800  J=1 ,130 
490  IARRAY( J)=0 
300  800  CONTINUE 
510  IF0UND=0 
520  RETURN 
530  END 
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10  INTEGER  I ARRAY ( 1 30  > 

20  CHARACTER*^  SN.HOLDSN 
30  CALL  INIT (1F0UND, IARRAY) 

40  READ(30, 10,END*500)H0LDSN 
50  10  F0RMAT(A15) 

60  50  READ(40, 100, END=499)SN, IDATE, I  DEMAND 
70  100  FORMAT (T1 ,A15,T67,I4,T79,I6) 

80  20  IF(SN.NE.HOLDSN)GO  TO  150 
90  IF0UND=1 

100  CALL  T ALL Y(SN, I DATE, IDEMAND, I ARRAY) 

110  GO  TO  50 

120  150  IFdFOUND.EQ.  1  )G0  TO  300 
130  IF(SN.GT.HOLDSN)GO  TO  200 
140  GO  TO  50 

150  200  URITE(6.250)H0LDSN 

160  250  FORMATC*  THIS  SN  NOT  ON  THIS  TAPE=\A15> 

170  GO  TO  450 

180  300  URITE(50,400)HOLDSN. ( I ARRAY ( I ) , I = 1 ,130) 

190  400  FORMAT (A1 5,5(72615)) 

200  CALL  INIT( IFOUND, IARRAY) 

210  450  REAH30, 1 0 , END=SOO ) HOLDSN 
220  GO  TO  20 

230  499  WRITE (6, 999) HOLDSN 

240  999  FORMATC  END  OF  TAPE=",A15) 

250  500  URITE(6,550) 

260  550  FORMATC  END  OF  SN  PROGRAM") 

270  STOP 
280  END 

290  SUBROUTINE  TALLY (SN, IDATE, IDEMAND, IARRAY) 

300  CHARACTER  SN*15 
310  INTEGER  IARRAY (130) 

320  IFdDATE.GE. 6276. AND. IDATE. LE.908?)G0  TO  600 
330  URITE(6,575)SN, IDATE, IDEMAND 

340  575  FORMATC  DATE  OUT  OF  BOUNDS*  ".AIS.SX.IA.SX.IS) 

350  RETURN 
360  600  CONTINUE 

370  IF( IDATE. GE. 6276. AND, IDATE. LE. 6366 )ISUB»(IDATE-6269)/7 
380  IF ( IDATE.6E.700t .AND.IDATE.LE.7365) ISUB*(IDATE-6994)77+13 
390  IF ( IDATE. GE. 8001. AND. IDATE. LE. 8365) ISUB*(IDATE-7993 >77+65 
400  IFdDATE.GE. 900 1 .AND. IDATE.LE. 9089) ISUB*(IDATE-8992>77+1  17 
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APPENDIX  D 

PROGRAM  TO  COMPUTE  FORECASTS  AND  PERFORM 
COMPARATIVE  ANALYSIS 
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U8*c 

1Z*=C 

138=C»*( 

STAR  PROCRAfl  CONFUTES  FORECASTS  US INC  NULTIPLE  TECHNIQUES 
AND  PERFORMS  BASIC  COMPARATIVE  ANALYSIS 

148* 

PROCCAR  ST AR  < I NPUT » OUTPUT » TAPE  1 • = I NPUT , Y  APE6=0UTPUT 1 

151* 

CQHHON  NSN (ZliNDATl 158) iIFSEI 158) iIFDEl 15(1 t 

168= 

1 IFAE ( 158) » IFHA4 (158) • IFNA8 ( 158 1 , 

17#= 

INS ( 1381 i INN ( 1581 tISEQ (35881 

181= 

READ  < 1 8 1 90 1 )  I 

181=981 

FORMAT (14) 

288= 

L=I 

218= 

DO  899  J=1,Z5 

228= 

READ (18.908) ISEQ ( I ) ,NSN ( 1 ) » NSN  (Z I » 

Z3I= 

1 (NOAT(IZ) ,12*1 , 1Z) » ISEQ ( I+l ! » (NDAT 

248= 

1 (IZ) .12*13.27) , ISEQ (1+2) . (NDATilZ). 

258= 

112=28.42), ISEQ (1+3). (NDAT IIZI. 

268= 

112=43, 57), ISEQ  (1+4MN0AT 112), 

Z7#= 

112=58, 72), 1SE0(I+5),(NDAT(IZI, 

288= 

112=73,87) ,ISEQ(I+6) , (NDAT 112) , 

Z9#= 

112=88, 182) ,  ISEQd+7) ,  (NDAT  ( 12) , 

3M= 

112=183, 117), ISEQ(I+8),(NDATII2), 

318= 

112=118,138) 

3ZI=9«8 

F0RHATU4, 11,  A18.A5, 1215,71/14, 15151, /I4.13I5I 

338= 

J5*L+8 

348= 

DO  185  J1*L,J5 

358= 

IFlISEQ(Jl).NE.Jl)  CQ  TO  898 

368=115 

CONTINUE 

378= 

CALL  SINCE! 

388= 

CALL  DOUBEI 

398= 

CALL  ADAPEX 

488= 

CALL  HQV4 

418= 

CALL  NOV8 

4Z8= 

CALL  MULTI 

438= 

CALL  ERROR 

448= 

L=L*9 

458* 

I*L 

468*899 

CONTINUE 

478= 

STOP 

488*898 

PRINT  975, ISEQ (Jl-l) 

498*975 

FORNAT (51 .Z8HCARD  OUT  OF  SEQUENCE, 41, 14! 

588* 

STOP 

518* 

END 
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53#*C  SUBROUTINE  SINCE!  COMPUTES  SINGLE  EIPONENTIAL  FORECAST 

54#=C  VALUES 

SSliClHHMMimHMWIMtHtilMIHHimiHtMmtitHHHHIH 

54#*  SUBROUTINE  SINCE! 

571*  COHNON  NSN (2) tNOAT < 158) * IFS£( 15# I > IFDE ( 15# ) * 

581*  1IFAE(150) »1FHA4(15#) *IFHA8(15#I » 

591*  INS (131) i INN (151) i ISEQ(3S«0 ) ■ 

4##=  IFSE(1)=MMTI1> 

411*  DO  2B8  1*1*129 

6Zi*  FSE*IFSE(1)+(.2»<NMT(I)-IFSE(I))) 

43#*  IFSE(I+1)*IFII(FSE) 

44#*Z##  CONTINUE 
45#*  RETURN 

44«*  END 

47#*C 

49f*C  SUBROUTINE  DOUSE!  CONFUTES  DOUBLE  EIPONENTIAL  FORECAST 

7##*C  VALUES 

71#*C 

72|:C*Hmt<liWtHHHH«<HHIiMHH««til«IH«litH«iillMIW 

73#*  SUBROUTINE  DOUBEI 

74#*  CONNON  NSN(Z) *NDAT(15#) *IFSE(15() tlFDE(lSI) * 

75#*  llFAEt  15# )  *  IFNA4U5#)  *  IFNA8  ( 15#)  * 

74f*  1NS(13#).INN(15#),ISEQ(35##) 

77#*  IFDE(l)*NOAT(l) 

78#*  IFOE(Z)--mTlZ) 

79#*  STC*(.f 

8##*  DO  3#f  1*2.129 

81#*  ATC=IFSE (I) -IFSE(I-l) 

82#*  STC=STC+ ( ,2» ( ATC-STC) ) 

83#*  IF (STC.LT .#. )  STC*#.# 

84#*C  5*1/. 2 

85#*  FDE* IFSE ( I ) + (5«STC) 

84#*  IFDE(I*1)*IFII(FDE) 

87#*3f#  CONTINUE 
88#*  RETURN 

89t*  END 
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921=0  SUBROUTINE  ADAPEI  COMPUTES  ADAPTIVE  EXPONENTIAL  FORECAST 

93#=C  VALUES 

940=C 

TShOMHMIHtHIlHHWtlHHHIIIIIIIIIMItUmmtlHiinttH 

941=  SUBROUTINE  ADAPEI 

971=  COMMON  NSN  (2) .NOAT ( 151) t IFSE 1 151 1 1 IFDE 1 151) i 
98f=  1IFAEI150) .1FHA4I150) .IFHA8I150) . 

99f=  INS (131) f inn ( 1S«) t ISEQ (3500T 

1»M=  IFAE(1)=NDAT(1) 

lflf=  ERT*0.0 

1020=  ENRT=«.I 

1130=  DO  400  1=1 *129 

1040=  IET=NDAT(I)-IFA£(I) 

1050=  IF(IET.E8.0.)  CO  TO  370 

1060=  CO  TO  380 

1070=370  AlPHA=0.0 

1080=  ERT=0.0 

1090=  ENRT=0.0 

1100=  CO  TO  390 

1110=380  ERT=(.2*IET)+(.8*ERT) 

1120=  IET=IABS(IET) 

1130=  ENRT=(.2*IET)+(.8*EHRT) 

1140=  ALPHA=ABSIERT/ENRT1 

1150=390  FAE= IFAE ( I ) ♦  I  ALPHA* I NDAT 1 1 1  - IFAE 1 1 ) 1 ) 

1160=  IFAE(I+1)=IFIX(FAE) 

1170=400  CONTINUE 

1180=  RETURN 

1190=  END 

1200=C 

1210=0************************************************************** 

12Z0=C  SUBROUTINE  NOV4  CONPUTES  4  TERR  RQVINC  AVERAGE 

1230=C  FORECAST  VALUES 

1240=C 

1Z50=C************************* **♦♦**+*♦* *************************** 
1260=  SUBROUTINE  N0V4 

1270=  COHHON  NSH (2). NDAT ( 150) r IFSE I 150) r IFDE! 151) r 

1280=  1IFAEU50)  ,IFHA4(150)  .IFHA8U50) . 

1290=  lHS(130Mffl(150)iISEQ<3500l 

1300=  DO  490  11=1.4 

1310=490  IFHA4(II)=NDAT!I1) 

1320=  IFHA4(1)=N0AT(4) 

1330=  DO  500  1=4.129 

1340=  FHA4=(NDATm*NDATU-l)*NDAT(J-Z)*NDAT<l-3)>'4 

1350*  IFNA4(I*1)=IFIX(FHA4) 

1360=500  CONTINUE 

1370=  RETURN 

1380=  END 
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1418=C  SUBROUTINE  ROV8  COMPUTES  8  TERM  MOVING  AVERAGE 

1*2»=C  FORECAST  VALUES 

1438=C 

1451=  SUBROUTINE  H0V8 

1461=  COMMON  NSN12) .NDATU58)  .IFSE1158)  .IFDEI158I > 

1478=  1 IFAE1158) » IFHA4 1 158 ) < IFMA8 ( 1581 . 

1488=  1HS1138) i INNI158) f 1SE013588) 

1498=  DO  598  11=1 .8 

1588=598  IFHA8 (Ill =NDAT (ID 

1518=  IFNA8(1)=N0AT(8) 

1528=  DO  688  1=8.129 

1538=  FNA8= (NDAT ( I I ♦NDAT ( I - 1 1 +NDAT ( I - 2 ) +NDAT ( I -3 1 + 

1548=  lNDAT(I-4)+NDAT(l-5)+N0AT(I-6HNDAT{l-7l)/B 

1558=  IFHA8(1+1)=IFI)UFNA8> 

1568=688  CONTINUE 

1578=  RETURN 

1588=  END 

1598=C»Mm#«#»m*M»Mm«Mmm*«m#n*mm#m«m«mm* 
1688=C  SUBROUTINE  MULTI  CONFUTES  MULTIPLE  MODEL 

1618=C  FORECAST  VALUES 

1628=C 

1648=  SUBROUTINE  MULTI 

1658=  COHNON  NSN(2)  .N0AT(1581.IFSEU58)  .IFDE1158I , 

1668=  11FAEU5I1  .IFRA4US8MFRA8U581 , 

1678=  lNS(138I.IHH<150).ISEfl<3588> 

1688=  DO  698  11=1.9 

1698=698  NS  (ID  =8 

1788=  00  728  1=9.129 

1718=  NVAR8=NDAT ( I ) +NDAT ( I- 1 ) +NOAT ( I-21+N0AT (1-31 

1728=  NVARl=(NVAR8-(IFSE(l)+IFSE(I-l)+IFSEtI-2)+ 

1738=  1IFSE4I-3) ) ) 

1748=  NVAR1=IABS(NVAR1) 

1758*  NVAR2=IA8$<NVAR8-(IFDEU)+IFDEII-l!+IFDE(I-2>+ 

1768=  lIFDE(I-3>>> 

1778=  NVAR3=IABS(NVA8B-lIFAEll)*IFAE(I-l>*IFAEtl-2D 

1788=  llFAEtI-3))) 

1798=  NVAR4-IABS (NVAR8* ( IFRA4< 1 1 +IFNA4 ( l- 1) ♦ IFNA4 ( I -2  > ♦ 

1888=  1IFNA4(I-3)D 

1818=  NVAR5-IABS (NVAR8- ( IFNA8 ( I 1+IFNA8 ( I - 1 ) ♦ 

1828=  lIFHA8(I-2)+IFAA8(l-3))l 
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83#* 

N*1 

184«* 

NVAR--NVAR1 

185#* 

IFINVAR2.LT.NVAR)  CO  TO  7(5 

186#*7#1 

IF (NVAR3.LT.NVAR)  CO  TO  796 

1871*7(2 

IF (NVAR4.LT. NVAR)  CO  TO  797 

188#*7#3 

IF (NVAR5.LT. NVAR)  CO  TO  718 

189#= 

CO  TO  719 

19##*7#5 

H=2 

191#* 

NVAR--NVAR2 

192#= 

CO  TO  791 

1930*706 

N*3 

19*#= 

NVAR=NVAR3 

195#* 

CO  TO  792 

1969=707 

N=* 

197#= 

NVAR=NVAR4 

198#* 

CO  TG  703 

199#=7#8 

H*5 

2999=71# 

HS(I)*N 

2919= 

IF(NSd).EQ.l) 

IHH(l+l)=IFSE(l*t) 

2929* 

1F(NS(I1.EQ.2) 

IHHll4l)*IFD£(l4l) 

2939* 

IF  (I9S  (I)  .EQ.31 

Iffl(I+l)*IFAE(l4l) 

29*9= 

IF (NS(I) .E8.4) 

IHN(I+1)*IFNA4(I+1) 

2959* 

IF(NS(I).EQ.5) 

IWl(l4l)*IFNA8(I+l) 

2069=729 

CONTINUE 

2070= 

RETURN 

2089= 

END 

ZI0#=C  SUBROUTINE  ERROR  MEASURES  THE  DECREE  OF  ERROR  AND  BIAS  IN 

2111*0  THE  FORECASTS 

2 120*C444ti44444444444444444#444444444444444444H444f  44444444  4444t##M 

213#=  SUBROUTINE  ERROR 

21*9*  COM  ON  NSN(2),NDAT(15#MFS£ll5»bIFDE(150t. 

2150*  1IFAE(15#)  »IFHA4(15#)»lFHA8d59) 1 

216#=  1H$!13#MMH(15#)»ISE8(35##) 

217#*  DIMENSION  ICTRU3#) 

218#*  SEN AD >9. 1 

219#*  DENAD*#.# 

229#*  AEMAD*#.# 

22  If*  A4NAD*#.# 

222#*  A8NAD*#.# 

223#*  SERE*#.# 

22*#*  HERE*#.# 

2251=  AER£*(.f 

226#*  A4RE*(.I 

227#*  A8R£*f.« 

228#*  SEPE*l.f 

229#*  HERE*#.# 

23##*  AEPE*«.f 

231#*  A4PE*#.# 

232#*  A8PE*f.l 

233#*  FNAD*#.# 

2341*  FRE*#.f 

235#*  FPE*#.# 
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3«*  DO  8#f  1*13.131 

2371*  ERR  *  NDAT(I)-IF$E(I) 

238#*  $EJM0*ABSIERRh$EHAD 
2398*  S£RE*ERR+SERE 
2W*  1F1ERR.EQ.8.)  CO  TO  757 

2<1«*  IFtNDAT (I) .EQ.C. )  CO  TO  751 
2«2f*  SEPE* (ABS (CRR/NOAT ( 1) ) *181, l+SEPE 
2431*  CO  TO  752 

2441*751  SEPE*  (ABS  (ERR/ 1 I  »1N .  l+SEPE 

2458*752  ERR*NDAT(1)-IFDE(I) 

244#1  OEIWD*ABSIERR)+OEIWD 

247#=  DERE=ERR+DERE 

248#=  IFIERR.EO.I.)  CO  TO  754 

2491=  IF(MDAT(I).EQ.|.)  CO  TO  753 

25M=  DEPE* (ABS (ERR/NDAT (1) 1*181, l+OEPE 

2511*  CO  TO  754 

252i*753  OEPE*  (ABS  (ERR/ 1 )  t  IN .  l+DEPE 

2531*754  ERR*NDAT(1)-IFAE(I) 

254#*  AENAD*ABS (ERR) +AEHAD 

2551*  AERE=ERR+AERE 

2541*  1F(ERR.EQ.I.)  CO  TO  754 
257#=  IF (NDAT (I).EQ.I.)  CO  TO  755 
258B*  AEPE* (ABS (ERR/NDAT ( I ) ) »1|| , ) +A£P£ 

25?#*  CO  TO  754 

Z4«*755  AEPE*  (ABS  (ERR/ 1 1  *1N .  1  +AEPE 

241#* 754  ERR*NDAT(I)-IFAA4(I) 

2428=  A4NA0*ABS(ERR1+A4IW1 

2431*  A4RE=ERR<44RE 

2448*  IF(ERR.E0.8.)  CO  TO  758 

2458=  IF (NDAT (I).E6.8.)  CO  TO  757 
2448*  A4PE* (ABS (ERR/NDAT (I) )*1#|. ) +A4PE 
2478=  CO  TO  758 
2488*757  A4PE*(ABS(ERR/1)*IN.)+A4PE 
2498*758  ERR*NDAT(IMFHA8(I) 

27N*  A8NAD* ABS  (ERR 1  +A8NAD 
2718*  A8RE*ERR+A8RE 
2728*  IF(ERR.E0.8,1  CO  TO  748 

2738*  IF(MAT(n.Ea.8.)  CO  TO  759 
2748*  A8PE* (ABS (ERR/MOAT (1)1*188.1 +A8PE 
2758*  CO  TO  748 

2748*759  A8PE*<A8S(ERR/1)»IN.)+A8PE 
2778*74#  ERR*NDAT(1)-I«mi) 

2708*  FIM0«AB$(ERR)+FMD 

2798*  FRE*ERR+FHE 

2888*  IF (ERR. EO. 8.)  CO  TO  745 

2818*  IF (NDAT (I). EO. 8)  CO  TO  741 

2828*  FPE* (ABS (ERR/NDAT (I) )*1##. ) +FPE 

2838*  CO  TO  745 

2848*741  FPE*(A8$(ERR/1)*1N,)+FPE 

2858*745  CONTINUE 

2848*88#  CONTINUE 
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•7#*  SEIMD=SE«AD/118. 

SERE=(SERE/118.»*t-l.#> 

28«--  SEPE=SEPE/118, 

DERAD=DERAD/U8. 

Z’t*=  DERE=(D£R£/ll8.)»(-j.0) 

mi=  0EPE»D£PE;il8. 

2930=  AERAD*AEAAD/H8. 

29<8=  A£RE=(AERE/118.)+(-l.0) 

2950=  AEP£=AEPEM18. 

2940=  A4«AD=A4M0/118. 

2970=  A4R£=!A4RE/lt8.)»(-lJI 

2980=  A4P£*A4PE/U8. 

2990=  A8HAD*A8!WD/118. 

300»=  A8R£=(A8R£/U8.I«(-1.|| 

3010*  A8PE=A8PE>  1 18. 

3«2i*  FRAD=FHAD)118. 

3831s  FRE=(FRE/118.)»(-1,0) 

3»4»s  FPEsFPE/118. 

3050=  ICTR(1)=1 

3848=  00  848  1=2,138 

3878=848  ICTR<I)=ICTR(I-1I+1 

3888=  WRITE (4,908)  NSAfl) ,NSN(Z1 

3898=900  F0RAATUH1,///5X,5HNSN:  ,2A10) 

3I#0=  13*1 

3110=  14=15 

3120=  ij.i 

3130=  DO  870  12=1,8 

3140=  IF1II.E8.4)  CO  TO  845 

3150=  CO  TO  841 

3140=845  URITE(4,912) 

3178=912  FORMAT ( 1H1 ) 

3180=  II=l 

3190=841  WRITE (4,901) (ICTR(I),I=I3,I4) 

3200=901  FORMAT  (//5I ,8H«»W£EH++ ,  4M3, 1 417) 
3210=  WRITE (4,902) (NDAT(I) ,1*13,14) 

3220=902  FORMAT </5l .4HACTUAL ,41 ,15,1417) 

3230=  WRITE(4, 903) (IFSE(I), 1=13,14) 

3240=903  F0RAAT(/5I,8H1-SIWCLE>2X,I5,14I7) 

3250*  WRITE(4, 904) (IFD£(I), 1=13,14) 

3240=904  FQRIWT(/5I,8W2>00UBLE, 21,15,1417) 

3270*  WRITE(4,905)(IFAE(I),I*I3,I4) 

3280*905  F0RMTt/5l,7H3-ADAPT,3M5, 1417) 

32?0*  WITE(4,904)(I«A4(I),I*I3,I4) 

3380=904  F0RlMT(/5J,7H4-4TERN,3J, 15,1417) 

3310=  WRITE(4»907)(IFNA8 (I), 1*13,14) 

3320*907  F0RMATC/5S,7H5-8TERM»3I,I5,14I7) 

3330=  WRITE(4, 908) (HS(I), 1*13,141 

3340*908  FORMAT ( /5I »4HNEn(CD ,41,15,1417) 

3350*  WRITE(4,909>  (Ilfflll) ,1*13,14) 

3340*909  F0RflAT(/5I,8HRULTIH00, 21,15, 1417) 

3370*  13*13+15 

3380*  14*14+15 

3390*  II*II+l 

3400*871  CONTINUE 
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♦II*  WRITE  (6*9»1HICT1M>!  1*13*131) 

3421*  WRITE (4»91f) 

3431*911  FORMAT  ( 1H-*  *9IX  *  3HHAD 1 8X  »4H8I AS 1 4X 1 4HRAPE ) 

3441*  WRITE(4»9iZMNOAT(I)»  1*13*  131) 

3451*  WRITE (4»9l3)(IFS£CI)tI*I3»13l) 

3441*  WRITE (6*911)  SEHAD » SERE » SEPE 

3471*911  F0RIWT(lH*i84Xi3Flt.41 

3481*  WR1TE(6»984)  CIFOEH)  *I-I3»13») 

3491*  WRITE (4i911)  DEHAD.DERE.OEPE 

35M*  WRITE(6.M5HIFfl£m.I*I3fl3»> 

3511*  WRITEI4r911)  AEHAD » AERE » AEPE 

3521=  WRITE(4*984)(IFHA4(I)it=13*13«> 

3531*  WRITE(4*911)  A4NAD*A4RE.A4PE' 

3541*  WRITE  (6t9i7HIFRA8(I)  1 1*13.1311 

3551*  WRITE(4.911)  AfiWADr A8RE * ASPE 

3541*  WRITE(4*988>  (NSU)*l*I3il3l) 

3571=  WRITE(6»9I9)  (UHI(t)  *1*13*131) 

3581*  WRITE14.911)  FWAB.FRE.FPE 

3591*  RETURN 

34H*  END 
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